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We developed 18 polymorphic simple sequence repeat (SSR) markers in pineapple (Ananas comosus) by us- 
ing genomic libraries enriched for GA and CA motifs. The markers were used to genotype 3 1 pineapple ac- 
cessions, including seven cultivars and 11 breeding lines from Okinawa Prefecture, 12 foreign accessions 
and one from a related species. These SSR loci were highly polymorphic: the 31 accessions contained three 
to seven alleles per locus, with an average of 4.1. The values of expected heterozygosity ranged from 0.09 
to 0.76, with an average of 0.52. All 31 accessions could be successfully differentiated by the 18 SSR mark- 
ers, with the exception of 'N67-10' and 'Hawaiian Smooth Cayenne'. A single combination of three markers 
TsuAC004, TsuACOlO and TsuAC041, was enough to distinguish all accessions with one exception. A pheno- 
gram based on the SSR genotypes did not show any distinct groups, but it suggested that pineapples bred in 
Japan are genetically diversed. We reconfirmed the parentage of 14 pineapple accessions by comparing the 
SSR alleles at 17 SSR loci in each accession and its reported parents. The obtained information will contrib- 
ute substantially to protecting plant breeders' rights. 

Key Words: Ananas comosus, genetic diversity, parentage, simple sequence repeat. 



Introduction 

The pineapple {Ananas comosus (L.) Merr.) is the most eco- 
nomically important edible plant of the family Bromeliaceae, 
which includes about 2,000 species, most epiphytic and 
many strikingly ornamental (Morton 1987). Pineapple is cul- 
tivated in most tropical and subtropical countries and in other 
regions with mild climates, ranking third in world produc- 
tion among tropical fruits, after banana and citrus (Botella 
and Smith 2008). Many pineapple cultivars are grown, dif- 
fering in characteristics such as plant and fruit size, flesh 
color and flavor and leaf margin type. Nearly all cultivars for 
commercial production are classified into a particular "type" 
category; examples include Cayenne, Queen, Maipure, Red 
Spanish, Singapore Spanish, Abacaxi and Cabezona (Wee 
and Thongtham 1991). 

In Japan, high-quality pineapple fruits can be produced 
only during the summer in the Ryukyu Islands, which stretch 
southwest from Kyushu to Taiwan (Republic of China), 
because the islands have a subtropical climate with mild 
winters and hot summers. Fruits harvested in winter are not 
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suitable for the fresh-fruit market, because of low tempera- 
tures during fruit maturation. Pineapple cultivation in the 
Ryukyu Islands and Okinawa Prefecture started in the 1920s 
or 1930s after immigrants from Taiwan brought pineapples 
to these islands and the pineapple canning industry was im- 
portant from the 1950s to the 1970s (Lin 1983, Watanabe 
1961). After the trade of processed pineapple fruits was lib- 
eralized in 1990, the proportion of pineapple production in- 
tended for the fresh-fruit market gradually increased. In 
2010, pineapple production in Okinawa was about 10,000 1 
(Ministry of Agriculture, Forestry and Fisheries Statistical 
Yearbook; http://www.maff.go.jp/j/tokei), 60% for fresh 
fruit consumption and the remainder for processing. 

Systematic pineapple breeding started in Okinawa Pre- 
fecture in 1989 and seven new elite cultivars have since been 
released for the fresh-fruit market. Several different types of 
cultivars (Cayenne, Queen, Maipure, Spanish and others), 
breeding lines and foreign accessions have been used as 
sources of specific characteristics {e.g., early ripening, high 
sugar content and low acidity). To establish effective breed- 
ing strategies, it is necessary to assess the genetic back- 
grounds of Japanese cultivars and other breeding materials 
by using molecular markers. In addition, it will be necessary 
to establish DNA profiling technique to protect new elite 
cultivars of pineapple. 
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In the plant variety protection (PVP) system in Japan, two 
main points were added by "Amendment of the Act in 2005" 
(http://www.hinsyu.maff.go.jp/en/about/overview.pdf) un- 
der the national policy for strengthening of intellectual prop- 
erty right. One was that coverage of plant breeders' rights 
(PBR) was expanded for products directly obtained from 
harvested material of the protected variety, because variety 
identification technique by DNA analysis has been devel- 
oped. Another was extension of the duration of PBR for 
30 years for fruit trees and woody plants. Therefore, DNA 
profiling technique would be important for protection of 
PBR in fruit species including pineapple. 

Up to now, several DNA profiling techniques have been 
used for cultivar identification and for evaluating genetic di- 
versity in pineapple. Ruas et al. (1995) used random ampli- 
fied polymorphic DNA (RAPD) markers to estimate the rela- 
tionships among four major pineapple cultivars in Brazil and 
Popluechai et al. (2007) analyzed nine pineapple cultivars in 
Thailand with 40 RAPD primers. Kato et al. (2004) used 
amplified fragment length polymorphism (AFLP) markers 
to evaluate 148 accessions of A. comosus and 14 of related 
species. Among the available molecular markers, simple se- 
quence repeat markers (SSRs, also known as microsatellites) 
provide a reliable method for cultivar identification because 
of their co-dominant inheritance and the abundance of al- 
leles per locus (Weber and May 1989). Wohrmann and 
Weising (2011) developed 18 EST-SSR loci that showed 
polymorphism in pineapple. However, DNA profiling in 
pineapple by SSR markers has been rarely studied. 

In this study, we developed 1 8 new genomic SSR mark- 
ers in pineapple by using an enriched genomic library ap- 
proach. We used them in cultivar identification, genetic 
diversity analysis, and parentage reconfirmation of 3 1 pine- 
apple accessions including Japanese cultivars, breeding lines 
and foreign accessions. 

Materials and Methods 

Plant materials and DNA extraction 

The 3 1 materials consisted of seven newly released cul- 
tivars from Okinawa Prefectural Agricultural Research 
Center Nago Branch (OPARC-Nago, Okinawa, Japan); 1 1 
breeding lines from OPARC-Nago; 12 foreign accessions 
introduced from the USA, Brazil, Taiwan and Australia and 
one related species, Ananas ananassoides (Table 1). All ma- 
terials were maintained at OPARC-Nago. Genomic DNA 
was isolated from young leaves by using a Genomic-tip 
20/G (Qiagen, Germany) as described by Yamamoto et al. 
(2006) or by using a DNeasy Plant Mini Kit (Qiagen, 
Germany) according to the manufacturer's instructions. 

SSR development 

We used genomic DNA of pineapple 'N67-10' to con- 
struct SSR-enriched genomic libraries for GA and CA mo- 
tifs by the method described by Nunome et al. (2006). The 
repeat-enriched genomic DNA was ligated into the pCR2.1- 



TOPO vector (TOPO TA Cloning Kit, Invitrogen, the 
Netherlands). Plasmid DNA was isolated and sequenced 
with the BigDye Terminator v3.1 Cycle Sequencing Kit 
(Applied Biosystems, USA) according to the manufacturer's 
instructions. Sequencing analysis was conducted with an 
ABI PRISM 3130x1 sequencer (Applied Biosystems, USA). 

A total of 384 plasmid sequences were obtained: 192 
from GA-enriched genomic libraries and 192 from CA- 
enriched libraries. Plasmid sequences with no insert were 
excluded from further analysis, and the minimum number of 
SSR repeats for marker development was set as eight repeats 
for di-nucleotide motifs of GA/CT or CA/GT. Primer pairs 
were designed with the Primer3 Web interface (Rozen and 
Skaletsky 2000, http://frodo.wi.mit.edu/primer3/input.htm). 
The general primer-picking conditions included a primer 
size of 20-25 bp (optimum 23 bp), a primer T m of 57-67°C 
(optimum 63°C), a maximum T m difference of 1°C, a primer 
GC content of 50%-60% (optimum 55%) and a product size 
range of 100-300 bp. 

SSR analysis 

SSR-PCR amplification was performed in a 10-|iL re- 
action mixture containing 5 |iL of GoTaq Master Mix in- 
cluding GoTaq DNA Polymerase (Promega, USA), 5 pmol 
each of forward primer (fluorescently labeled with Fam, 
Vic, or Ned) and reverse primer (unlabeled), and 5 ng of ge- 
nomic DNA. DNA was amplified in 35 cycles of 94°C for 
1 min, 55°C for 1 min and 72°C for 2 min and a final exten- 
sion of 10 min at 72°C. The nucleotide sequence of "gtttctt" 
was added to the 5' end of reverse primers as pig-tailing 
(Brownstein et al. 1996), in order to promote adenylation 
and facilitate accurate genotyping. The amplified PCR prod- 
ucts were separated and detected in a PRISM 3100 DNA 
sequencer (Applied Biosystems, USA). The sizes of the am- 
plified bands were scored against internal-standard DNA 
(400HD-ROX, Applied Biosystems, USA) by GeneScan 
software (Applied Biosystems, USA). 

Data analysis 

Using the CERVUS v. 2.0 software (Marshall et al. 1998) 
and MarkerToolKit v. 1.0 (Fujii et al. 2008), we estimated 
the expected heterozygosity (H E ) at single-locus SSR mark- 
ers in the tested pineapple cultivars. H E was calculated using 
an unbiased formula from allele frequencies as 1 - Sp/ 2 
(\<i<m), where m is the number of alleles at the target lo- 
cus and pi is the allele frequency of the rth allele at the target 
locus. 

Parent-offspring relationships were tested by comparing 
the SSR alleles in each accession with those of its reported 
parents; the data were analyzed using the MARCO software 
(Fujii et al. unpublished). MinimalMarker software (Fujii et 
al. 2007) was used to identify minimal marker subsets to dis- 
tinguish all cultivars and to find identical genotypes generat- 
ed from the 18 SSR markers for the 31 accessions. 

A phenogram of the 3 1 accessions was constructed by us- 
ing the unweighted pair-group method with using arithmetic 
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mean (UPGMA) based on the similarities between geno- 
types estimated by Dice's coefficient: D c = 2n xy l{n x + n y ), 
where n x and n y represent the number of putative SSR alleles 
for materials X and Y, respectively, and n v represents the 
number of putative SSR alleles shared between X and Y. 
NTSYS-pc v. 2.1 software (Rohlf 1998) was used to visual- 
ize the phenogram. 

Results 

SSR marker development 

We sequenced 384 plasmid clones from GA- and CA- 
enriched genomic libraries of 'N67-10' (192 clones from 
each library). After exclusion of clones with no inserts, am- 
biguous nucleotide sequences, no repeat motifs, and duplica- 
tion, 110 sequences that contained at least eight repeats of a 
di-nucleotide motif remained from the GA-enriched geno- 
mic library. These sequences contained 8 to 49 repeat motifs 
of (GA)/(CT), with 16.8 on average. The average insert size 
of the clones obtained was about 202 bp, ranging from 84 to 
443 bp. Ninety-eight sequences were obtained from the CA- 
enriched library containing 8 to 31 repeats of (CA)/(GT), 
with an average of 12.8. The average insert size of the clones 
obtained was about 255 bp, ranging from 122 to 552 bp. A 
total of 42 primer pairs were designed with the Primer3 pro- 
gram: 23 for GA repeats and 19 for CA repeats. 

We screened and evaluated the 42 SSR primer pairs by 
using eight pineapple accessions: 'N67-10', 'Cream Pine- 
apple', 'Julio Star', 'Summer Gold', 'Yugafu', Okinawa 
No. 17, 'Soft Touch' and Okinawa No. 19. Out of the 42 
SSR marker candidates, 24 were excluded from further anal- 
ysis because of no amplification or unstable amplification of 
the target band. The remaining 18 markers (7 with GA- 
repeat motifs and 1 1 with CA-repeat motifs) were used for 
SSR analysis of all 31 pineapple accessions (Table 2). 
Among the seven GA-repeat SSR markers, five showed per- 
fect repeats of a GA motif, whereas TsuACOlO and 
TsuAC013 had combined motifs of (GT) 14 A(AG) 12 and 
(AGAGAT) 3 (AG) 12, respectively. Out of the 11 CA-repeat 
SSRs, nine showed perfect repeats of a CA motif, whereas 
TsuAC018 and TsuAC023 had an interrupted CA motif of 
(CA)i 0 A(AC) 9 and a combined motif of (CA) 10 (TA) U , re- 
spectively. 

Genetic identification of pineapple 

We identified 74 putative alleles in the 3 1 pineapple ac- 
cessions with the 18 SSR markers (Table 2). The number of 
alleles per locus ranged from three at five of the loci 
(TsuAC024, TsuAC028, TsuAC035, TsuAC038 and 
TsuAC039) to seven at TsuACOlO, with an average value of 
4. 1 (Table 2). The expected heterozygosity (H E ) ranged from 
0.09 at TsuAC019 and TsuAC026 to 0.76 at TsuACOlO, 
with an average value of 0.52. 

The 3 1 pineapple accessions could be successfully differ- 
entiated from one another by the 18 SSR markers, with the 
exception of 'N67-10' and 'Hawaiian Smooth Cayenne' 



(Table 3 and Fig. 1). A single combination of three markers 
(TsuAC004, TsuACOlO and TsuAC041) was enough to 
distinguish 30 accessions (all except for 'N67-10' and 
'Hawaiian Smooth Cayenne') on the basis of at least one dif- 
ference in SSR genotype. Furthermore, ten marker subsets 
consisting of six SSR markers each (e.g., TsuAC004, 
TsuAC008, TsuACOlO, TsuAC030, TsuAC039 and 
TsuAC041) could differentiate 30 accessions on the basis of 
two or more SSR genotype differences. 

We constructed a phenogram of the 3 1 accessions based 
on SSR analysis (Fig. 1). The accession belonging to the re- 
lated species A. ananasoides was clearly separated from the 
other 30 pineapple accessions. These 30 pineapple acces- 
sions were not separated into distinct groups but seemed to 
be mingled together, and there was little relationship be- 
tween the cultivar types and the genetic distances based on 
the SSR analysis. Although some accessions were used repre- 
senting several types, i.e., 'Cream Pineapple' (Maipure type), 
'McGregor ST-F (Queen type), 'Red Spanish' (Spanish 
type), 'Seijyo Cayenne' and 'Hawaiian Smooth Cayenne' 
(Cayenne type), no distinct groups were found. 

Parentage analysis 

We examined parent-offspring relationships of 1 5 pine- 
apple cultivars and breeding lines bred by OPARC-Nago by 
using 17 SSR genotypes of the 18 loci (Table 3, all except 
for TsuAC007). SSR analysis suggested that 'N67-10' had 
been selected and bred from a sport or mutant of 'Hawaiian 
Smooth Cayenne' (Ikemiya et al. 1984). The parentage anal- 
ysis was conducted by comparing SSR alleles in the acces- 
sions and their reported parents. The parentage of ten acces- 
sions ('Soft Touch', 'Summer Gold', 'Yugafu', 'Julio Star' 
and Okinawa No. 9, 13, 17, 19, 22 and 24) was reconfirmed: 
in each of these accessions, the SSR alleles at each locus 
were consistent with one allele being contributed by each of 
the reported parents. Parent-offspring relationships were 
also reconfirmed between Okinawa No. 21 and 'Yugafu' 
and Okinawa No. 23 and 'Julio Star'; but in each case, the 
second parent was not available for testing. A discrepancy at 
one SSR locus TsuAC018 was found for 'Gold Barrel' and 
its reported parents 'Cream Pineapple' and 'McGregor ST- 
1', which might have been caused by a mutation, otherwise 
the existence of a null allele, i.e., 109/null, 109/121 and 121/ 
null genotypes for 'Gold Barrel', 'Cream Pineapple' and 
'McGregor ST-F, respectively. In the case of Okinawa 
No. 20, the SSR data were inconsistent with the reported par- 
entage ('Yugafu' x 'Summer Gold') but suggested another 
possible set of parents ('Yugafu' x 'N67-10'). 

Discussion 

Interest in plant breeders' rights is increasing worldwide. 
The International Union for the Protection of New Varieties 
of Plant (UPOV; http://www.upov.int/index_en.html) pro- 
motes the development and use of effective systems of plant 
variety protection. The BMT (Biochemical and Molecular 
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Table 2. Characteristics of 18 nelwly developed SSR markers in pineapple 



SSR locus 
Accession nos." 


Primer sequence (5 '-3') 


Repeat motif 


Target size 
(bp) 


No. of 
allels 


Heterozygosity 


TsuAC004 


F: ATGTTGGTCAAAGGGCTGTT 


(AG) 16 


144 


5 


0.67 


AB7 16708 


R: gtttcttTCATGATCACACTGGAGATTTG 










TsuAC007 


F: GCAGCGGTAAGATCTGCTTT 


(GA)21 


102 


4 


_4 


AB716709 


R: gtttcttTCCTTCTCTCCACCTCTTCATT 










TsuAC008 


F: GAAATGGTACTGCTTCACTGTTC 


(GA)16 


173 


5 


0.71 


AB716710 


R: gtttcttATACGGGGAAATAGGCACAA 










TsuACOlO 


F: TGAGTTGTGTCATTGTGTGTCA 


(GT)14A(AG)12 


207 


7 


0.76 


AB716711 


R: gtttcttGGGGGTCTCCATACATTTTT 










TsuACOB 


F: TTATGCAGGAAAATAGGGGG 


(AGAGAT)3 (AG) 1 2 


139 


4 


0.55 


AB716712 


R: gtttcttC ATGC ATC ATA AATTC GTGTCC 










TsuAC018 


F: GCATCGATCTCCATGCAAAC 


(CA)10A(AC)9 


120 


5 


0.59 


AB716713 


R: gtttcttAAAGGAAACAAGGAGGATGTGA 










TsuAC019 


F: TTCATCCTATGGTTTCCCCA 


(AC) 13 


177 


4 


0.09 


AB716714 


R: gtttcttGTGGGTTCAACTGAGTAGCAAT 










TsuAC021 


F: AATCAAAGTGATTCCCCTTCC 


(CA)21 


141 


4 


0.50 


AB716715 


R: gtttcttTCTGACATAGGGCTTGCACA 










TsuAC023 


F: TCGAAAAGAGGATGCTGGAT 


(CA)10(TA)11 


143 


5 


0.73 


AB716716 


R: gtttcttTCCGCAGTGTAGGCATGTAA 










TsuAC024 


F: GTCGCCAATCAAATTCCAGT 


(AC)9 


126 


3 


0.52 


AB716717 


R: gtttcttCTCACGAAACATGAATCACCA 










TsuAC026 


F: GGGATTAACTTTTCCAGGGG 


(AC)8 


200 


4 


0.09 


AB716718 


R: gtttcttTTGGATTCCTCGTTTGCATT 










TsuAC028 


F: TGACACCATAGAGGAGGGGT 


(AC)8 


220 


3 


0.57 


AB716719 


R: gtttcttGCTCAAGGACAATCCACCAT 










TsuAC030 


F: GAGAGAGAAAAGAGTTTCGACAG 


(AG)27 


149 


4 


0.43 


AB7 16720 


R: gtttcttCTTCAAAATGGTCTAACGTACC 










TsuAC035 


F: TTCCTAGCCAACACTACTACAGA 


(GA)9 


96 


3 


0.45 


AB7 16721 


R: gtttcttTGCAGCTTCTTTTCCTGGTT 










TsuAC038 


F: TTGCAGCAAACCAAGTCAT 


(AC) 11 


327 


3 


0.48 


AB716722 


R: gtttcttGGAGGTGTAGTCAATTAGGAGAA 










TsuAC039 


F: CCCTGTATGGGTAGCATTGAA 


(AC)8 


91 


3 


0.54 


AB7 16723 


R: gtttcttAAAAGGTATCACGAAAGCGA 










TsuAC040 


F: AAATTCTCTTCATGCACACG 


(AC)8 


99 


4 


0.61 


AB7 16724 


R: gtttcttTGCTTCATGAGATCTAAACTGG 










TsuAC041 


F: CTCTCTTATGGCACAACCCTG 


(AC) 11 


279 


4 


0.58 


AB7 16725 


R: gtttcttCCTGGTGAGTAATCTATATGCTG 










Average 








4.1 


0.52 



a DDBJ accession numbers 

b Heterozygosity not evaluated because of the existence of null allele 



Techniques, and DNA-Profiling in Particular) working 
group of UPOV agreed to establish DNA profiling tech- 
niques for protecting plant breeders' rights. In Japan, ten 
DNA profiling manuals for major crops, including rice 
(Oryza sativa L.), kidney bean (Phaseolus vulgaris L.), 
adzuki-bean (Vigna angularis (Willd.) Ohwi & Ohashi.), 
strawberry (Fragaria x ananassa Duchesne), sweet cherry 
(Primus avium L.) and Japanese pear (Pyrus pyrifolia 
Nakai), have been released on the Plant Variety Protection 



home page (http://www.hinsyu.maff.go.jp/en/en_top.html). 
Among them, methods for the use of SSR markers in fruit 
species are given for sweet cherry and Japanese pear. 

Genome-specific CAPS markers had greatly contributed 
to prevent an infringement of breeder's rights in strawberry 
(Kunihisa 2011). Packed strawberry fruits imported from 
Korea labeled as 'Nyoho' were identified to be a mix of 
'Redpearl' and 'Sachinoka' by using CAPS markers 
(Kunihisa et al. 2005). After admonition to dealers and 
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Fig. 1. Phenogram of the 31 pineapple accessions evaluated in this study. The phenogram was produced using the UPGMA method based on 
Dice's coefficient. 



action to the court, the volume of illegally imported straw- 
berry fruits sharply decreased. On the case of sweet cherry, 
'Benishuho' in which Yamagata Prefecture holds the breed- 
er's right, was unlawfully taken out overseas by an Australian 
citizen residing in Tasmania, producing and selling fruits. 
Thus, Yamagata Prefecture established DNA profiling 
system for fruit tissues using SSR markers and lodged a 
criminal complaint against the exporters (Takashina et al. 
2007, 2008). The other recent breeder's rights infringement 
cases, DNA profiling methods were developed for rush, 
kidney bean, adzuki-bean and etc. In this study, the SSR- 
based identification system for major pineapple cultivars in 
Japan was developed, which will contribute greatly to pro- 
tecting plant breeders' rights. 

Kato et al. (2004) reported that major cultivar types such 
as Cayenne, Spanish and Queen, could not be distinctively 
separated based on AFLP analysis using 148 accessions of 
A. comosus and 14 of related species. In this study, the 31 
pineapple accessions were not clustered into distinct groups 
in a phenogram constructed from the results of the SSR anal- 
ysis. Therefore, it is considered that cultivar types have been 
classified based on morphological similarity, and that DNA 
analysis was not in good accordance with morphological 
classification. Further DNA analysis will help us establish 
an accurate classification system. Our results suggested that 
abundant genetic variation existed within cultivars and 
breeding lines in Japan and foreign accessions. Discrete 
DNA profiling of pineapples by SSR markers will be uti- 
lized for cultivar protection systems. 

Accurate information about parent-offspring relation- 
ships is necessary for efficient breeding programs. SSR 
markers have been used for parentage analyses of grapes 
(Bowers and Meredith 1997, Bowers et al. 1999, Sefc et al. 



1997), peaches (Testolin et al. 2000, Yamamoto et al. 2003), 
apples (Kitahara et al. 2005, Moriya et al. 2011) and pears 
(Sawamura et al. 2008). We examined parent-offspring rela- 
tionships of 15 pineapple cultivars by using 17 SSR loci and 
reconfirmed the parentage for 14 of the cultivars. In the case 
of Okinawa No. 20, a candidate for its true parent combina- 
tion was identified. 
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